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ABSTRACT 

The composition of the nuclear component of the cosmic radiation 
has been under study for some time, because it can yield information 
concerning the source, propagation, and confinement of cosmic rays 
within the galaxy. We present here the first comprehensive satellite 
measurement of cosmic-ray composition and spectra in the charge range 
2 < Z s 10 using the geomagnetic field as a rigidity analyzer through 
the entire range of vertical cutoffs. The results indicate that the 
spectra of all the elements in the observed range are similar, and thus 
that various ratios of elemental abundances are nearly independent of 
rigidity over the range 2 s P s 15 GV. Calculations of the propagation 
of cosmic rays through the interstellar and interplanetary media pre- 
dict that there should be a variation with rigidity of ratios of 
various elements, because of the charge-dependent effects of ioniza- 
tion of the interstellar gas by the cosmic rays. The absence of this 
variation can be explained by assuming a rigidity-dependent confinement 
of the cosmic rays in the galaxy. From the ratios Be/0 and B/0 we 

obtain a leakage mean free path decreasing from -vlO g/cm at ^2 GV to 
2 

^<6-7 g/cm at 10-15 GV, assuming that the interstellar gas is composed 
of hydrogen and helium in a 10:1 ratio. Ratios of other cosmic-ray 
elements are consistent with such a decreasing path length. The con- 
finement time of the cosmic rays in the galaxy can be estimated from 
the ratio Be/B, because of the radioactive decay of ^°Be to 10 B. The 
data indicate a confinement time of -|o®‘ 5±0 ' 5 years, which may be 
rigidity-dependent. 
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I. Introduction 

Ever since the discovery ( Freier et al. , 1948) of nuclei heavier 
than protons in the cosmic radiation, considerable effort has been 
expended in making measurements of the elemental composition and the 
energy spectra of this component of the cosmic rays, and in developing 
models to account for the observations. The current picture of the pro- 
duction of the heavy nuclei can be outlined as follows. A source, or a 
number of sources, injects into the galactic medium ,a flux of high 
energy nuclei. The composition of this flux depends on the composition 
of the source region, but it is presumed to contain primarily protons, 
a-particles, and those nuclei (such as 12 C, 16 0, 20 Ne, 24 Mg, 28 Si, 56 Fe) 
which are expected to be produced abundantly in various thermonuclear 
processes which presumably have produced the source material. These high- 
energy particles propagate within the galaxy, guided by the interstellar 
magnetic field. Because of random irregularities in this field, the par- 
ticles diffuse throughout the galaxy, and may escape when they reach the 
boundary. During their travel, they may interact with atoms of the 
interstellar gas, which is composed mainly of hydrogen. They may ionize 
the atoms of the gas, in which case they lose energy to the atomic elec- 
trons, and they may suffer nuclear interactions in which they are trans- 
muted into lighter nuclei. By this latter process, nuclei are formed 
which may not have been present at the source. Thus, by measuring the 
composition near earth, we can deduce some of the characteristics of the 
propagation, and, given enough data and adequate models, we can also 
learn something about the sources. For instance, a knowledge of the con- 
finement time and the energy density of the cosmic rays provides 
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information on the overall source strength. 

Recently, the advent of many new technologies has brought this field 

to the point where detailed comparisons of theory and experiment are 

possible. Some of the recent developments which have been very useful 

are: solid-state detectors (silicon and large-area crystal and plastic 

scintillators), rugged, sensitive photomultiplier tubes, improved 
v 

detectors for Cerenkov radiation, and superconducting magnets, as well 
as large, high-altitude balloons and, of course, artificial satellites. 

An excellent review of the recent status of the measurements and theoret- 
ical implications was published by Shapiro and Silberberg (1970). 

Most of the measurements of composition and spectra have been made 
by experiments flown on high-altitude balloons ( e.g ., Webber and Ormes, 
1967; Cory don-Peter sen et at. , 1970; Smith et at., 1973; and many others). 
Experiments of this type can be quite large and massive, providing the 
large collection factors needed to obtain good statistics from the rather 
low particle fluxes, and allowing determination of particle energies or 
rigidities up to quite high values [e.g., rigidity ;> 100 GV). However, 

these experiments are shielded from direct exposure to the cosmic rays 
2 

by several g/cm of overlying atmosphere. Ionization and nuclear inter- 
actions in the air above the experiment produce distortions in the 
spectra and composition of the radiation, requiring corrections which may 
introduce a certain amount of uncertainty into the results. 

This problem can be avoided by making the observations on satellites 
outside the atmosphere. The difficulties here are imposed by the neces- 
sarily small size of most satellite experiments. In order to obtain 
adequate counting statistics, data must be taken for periods on the order 
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of months or years, requiring considerable stability and reliability of 
the instrument. Also, it is difficult to make determinations of energy 
above 1 GeV/nucleon, due to the small size and mass available for 
experiments. As a result, most satellite measurements are restricted 
to energies below a few hundred MeV/nucleon ( e.g ., Balasubrahnanyan 
et at. , 1966; Comstock et at., 1966). This can be extended to ^ 1 
GeV/nucleon by using a Cerenkov detector in addition to the usual dE/dx 
and total- E detectors { Leaniak and Webber, 1969; Teegarden et at., 1970; 
Mason , 1972). 

It has been known for quite some time ( Orear et at., 1949; Stormer, 
1955 and references therein) that the geomagnetic field prevents some 
particles from reaching the top of the atmosphere. The simplest model 
of the process, using a dipole field and ignoring the presence of the 
solid earth, predicts that particles having a magnetic rigidity 
(momentum/charge) below a certain cutoff value cannot be incident at the 
surface of the earth. For vertically incident particles, the cutoff is 
given simply by 

P c = PqCOs 4 x (1-1) 

where X is the magnetic dipole latitude, and P Q = 15 GV (rigidity is 
commonly defined as pc/Ze and measured in volts; 1 GV = 10^ volts). 

Thus the geomagnetic field can be used to determine threshold rigidities 
up to 'v 15 GV, depending on latitude, enabling one to determine spectra 
up to that rigidity. The actual calculation of cutoff rigidities is 
considerably more complex, however, than is indicated by Equation 1-1 . 
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The method used is described in Appendix A. 

We describe in this thesis the first satellite measurement of cosmic 
ray elemental composition using the geomagnetic field as a spectral 
analyzer, through the entire range of vertical cutoffs from zero to 
^15 GV. Techniques of data analysis are described which yield a charge 
resolution of ^0.2 charge unit. Integral rigidity spectra are presented, 
as well as ratios of various elements. The results are discussed in 
relation to a steady-state model for propagation and confinement of 
cosmic rays within the galaxy. 
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II. Experiment 

A. Spacecraft 

The Caltech Solar and Galactic Cosmic Ray Experiment ( Althouse 
et oil., 1967) was launched into a polar orbit aboard the 0G0-6 spacecraft 
on June 5, 1969. The initial orbit parameters were: perigee 397 km, 

apogee 1098 km, inclination 82°, period 99.8 min (0G0 Experiment 
Bulletin F-143). The data used in this thesis were taken during the 
period June 7, 1969 to August 28, 1970. 

The low-altitude polar orbit of 0G0-6 provided extensive sampling 
of the entire range of vertical geomagnetic cutoff rigidities up to 
M5 GV. In addition, the spacecraft was oriented so the Caltech experi- 
ment is always facing radially away from the earth. 

B. Cerenkov Telescope 

At moderately high energies, i.e., above a few hundred MeV/nucleon, 
a useful technique for charged-particle measurements is a simultaneous 
determination of ionization energy loss in a thin detector and of 
Cerenkov radiation in a material of suitable index of refraction, n. 

If the particle velocity 3 is greater than 1/n, then the number of 
Cerenkov photons emitted is proportional to Z [l-($n) ], where Z is 
the charge. The particle's energy loss in a thin detector is approxi- 
mately proportional to Z /3 . Thus a simultaneous measurement of these 
two quantities can be used to determine both the charge and velocity of 
the particle. 
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A cross section of the Cerenkov Telescope is shown in Figure II-l. 
D1 ' and D2‘ are totally depleted Au-Si surface barrier detectors, 
nominally 1 mm thick, with diameters of ^2.34 cm. D3 1 is the 1 cm thick 
fused silica window of a 2-inch photomultiplier tube. D4' is a plastic 
scintillator anticoincidence shield. A D2 , D3'D4' coincidence is 
required for event analysis. Particle energy loss (-dE/dx) is measured 
in D1‘ and D2' , and pulse-height analyzed over a range from % 100 keV to 
^0 MeV. This covers the expected energy-loss range from minimum- 
ionizing, singly charged particles to just above the value for minimum- 
ionizing magnesium nuclei. The fused silica of D3‘ has an index of 
refraction of about 1.6 in the near ultraviolet, dropping to about 1.45 

■ V 

in the visible. This corresponds to a Cerenkov threshold of about 
300 MeV/nucleon. 

The calculated response of the detectors to particles of various 
charges and energies is shown in Figure 1 1-2 . Here energy loss is cal- 
culated from {Evans, 1955) 

4 2 

47re Z NZ o 99 

dE/dx ^ [In (2m Jr/ 1) - ln(l - p ) - fT] (II-l) 

mj* e 

e 

where Z is the charge of the incident particle, V = pc' is its velocity, 

1 is the charge of the absorbing material (Z a = 14 for silicon), m is 

a 3 c 

3 

the mass of the electron, N is the number of absorber atoms per cm , 
and I is an empirical "excitation potential" whose value for silicon is 
M70 eV. For silicon, this can be rewritten 

i 

dE/dx = (.359 MeV/cm) (Z 2 /g 2 ) [8.70 + ln(e 2 y 2 ) - B 2 ] (11-2) 
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Figure I 1-1 

Cross-sectional view of 060-6 Cerenkov Telescope. 
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Figure II-2 

V 

Predicted response of 060-6 Cerenkov Telescope, 
light output is normalized so a singly charged, 
relativistic particle gives unit light output, 
using a value of 1.58 for index of refraction. 


Cerenkov 
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RELATIVE CERENKOV OUTPUT 
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where y is the Lorentz factor. 

Cerenkov output is calculated from [Evans, 1955) 

d 2 N/dxdv = 2ttciZ 2 c‘ 1 / 1 - -*4 ) (II-3) 

l en 2 (v) / 

2 

where d N/dxdv is the number of photons emitted per unit path length in 
the frequency interval v to v + dv, a =1/137, and n(v) is the index of 
refraction of the material. If we include the response of the photo- 
multiplier tube as e(v) = number of photoelectrons/photon, then 

K = 27raZ 2 C _1 AX J [1 - B‘V 2 (v)] e(v)dv (II-4) 

where the integration is done over all frequencies for which the inte- 
grand is greater than zero, and K is the total number of photoelectrons 
emitted from the photocathode. A reasonably good approximation for the 
0G0-6 photomultiplier is e =.18 for vacuum wavelengths from ^170 nm to 
^520 nm and zero elsewhere, and n = 1.58 in this range. Using 
Ax = 1 cm, we get 

K = 345 (1 - .40/b 2 )Z 2 photoelectrons ( I 1-5 ) 

for normally incident particles of charge Z. For fully relativistic 

p 

(g = 1), downward-moving particles, this is K = 200Z photoelectrons. 
The face of the silica window is painted black to discriminate against 
upward-moving particles ( Wenzel , 1968). 
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C. Calibrations 

Three types of calibrations are available to assist in the inter- 
pretation of data from this experiment. Electronic calibrations on the 
ground provide detailed information on pulse-height analyzer (PHA) 
linearity, temperature sensitivity, and stability. Calibrations using 
artificially accelerated protons (Caltech Tandem Van de Graaff) and 
electrons (Caltech Synchrotron), and natural muons, alphas, and 8- 
particles yield information on detector diameters, thicknesses, and 
uniformities. These calibrations are described in detail elsewhere 
( Murray 3 1971). The electronic calibrations are summarized here to 
illustrate the correspondence between detector outputs and PHA channels. 
Figure I 1-3 shows the 25 °C calibrations for DV and D2‘, and Figure I 1-4 
shows the same for D3‘. 

This experiment was designed with two identical dE/dx detectors in 
order to minimize the effects of statistical fluctuations in the energy- 
loss process ( Symon, 1948), and to reduce background from a number of 
sources. However, the overall responses of the two detectors (including 
amplifiers and PHA's) are not exactly identical, as can be seen from the 
calibrations shown in Figure II-3. These calibrations were confirmed by 
using the relativistic particle peaks for DV and D2', in a manner 
similar to that indicated below. In order to compare the two detectors, 
a correction was made to the DV response by defining 


Dl" = DV + .075 (DV - 190) 


( I 1-6 ) 
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Figure II-3 

Laboratory calibration of DV and D2' at 25°C. Piecewise 
linear transfer characteristic, explained in Althouse 
et at. (1967), is used to achieve large dynamic range 
without sacrificing resolution at low signal levels. 
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Figure II-4 

Laboratory calibration of D3‘ at 25°C. PHU (Pulse Height 
Unit) is a calibrating unit corresponding to a photomul- 
tiplier output of about 25 picocoulombs. One pi co- 
coulomb corresponds to a relative light output of about 
unity in Figure I I -2 . 
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for DV > 127. For DT = 125 to 127, Dl" = 125. The validity of this 
correction was checked by examining the means of the Dl" - D2' distri- 
butions as discussed in Chapter III. 

The most important calibration is provided by actual flight data 
accumulated at locations where the vertical geomagnetic cutoff rigidity 
(P ) was sufficiently high 6.4 GV) to ensure that only particles 

v 

having 3 = 1 entered the detector telescope. Note from equations ( 1 1-2 ) 
and ( 1 1-5 ) that, for a fixed value of 3, both energy-loss and Cerenkov 
outputs are linear in Z . Therefore, we expect particles having 3 - 1 

v 2 

to produce data points clustered along a line with aE“C«Z . Figure I 1-5 
illustrates the distribution in (Dl" + D2')/2 versus D3 ' space of 3 - 1 
events selected according to consistency ‘criteria discussed in a later 
section. Since the individual charge peaks are clearly separated, it 
is a simple matter to count up the number of events of a given. charge, 
as a function of either AE channel or l channel. While this procedure 
is subject to some error in charge assignment for events far from the 
peaks, the error is not large, and in any case is not important for the 
present purpose of determining the locations of the peaks. A more 
sophisticated charge determination technique is discussed in Chapter III. 

The results of this procedure for the boron and nitrogen peaks are 
shown in Figure I 1-6 . Similar plots were made for all the peaks, and 
the results for 3<Z<10 are shown in Figure I 1-7 . As expected, we have 
aE<*Z and CccZ . The lines fitted to the peak positions form the basic 

i 

calibration of the instrument. This calibration may be expected to be 
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Figure II-5 

Scatter plot of high-rigidity events showing positions 
of peaks for relativistic particles. Sample includes 
only events with P > 6.4 GV to ensure that only par- 

V 

tides having 3 = 1 contribute to the distributions. 
Vertical axi£ represents average of two energy-loss 
measurements; horizontal axis is Cerenkov output. 
Numbers on plot represent number of events. Letter 
code is logarithmic, with each letter representing 
a factor of 10 0#1 , i.e., A - 10, F - 32, K - 100, etc. 
Events selected as discussed in Chapter III. Light 
lines delimit preliminary charge assignment used only 
for determination of locations of peaks. 
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Figure II-6 

Events selected from boron and nitrogen distributions 
of Figure I 1-5. Accurate charge identification was 
not attempted for these plots, the purpose of which 
is to locate the centers of the distributions. 
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Figure I 1-6 
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Figure I 1-7 

Peak positions of distribution shown in Figure II-5, 
obtained from plots such as those in Figure II-6. 
Uncertainties in locations of peaks are sO.2 channel 
for Z < 8. Poor statistics make locations more 
uncertain for Z > 8. 
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temperature sensitive, due to variations of the characteristics of the 
electronics with temperature. However, the temperature of the space- 
craft -Z door, on which the experiment is mounted, was held at 
21.3 ± 1.9°C by the spacecraft Thermal Control System. Comparison of 
data taken when the temperature was s2Q°C and ;>23°C indicate a shift of 
£1/2 channel in the position of the oxygen peak (Figure I 1-7), so all 
data discussed in this thesis were combined without regard to 
temperature. 
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III. Data 

A. Event Selection Criteria 

In order to take advantage of the redundant energy-loss measure- 
ments provided by Dl' and D2 1 , we need to understand the details of the 
overall response functions of the detector systems. The energy-loss 
distributions in Dl 1 and D2 1 are highly asymmetrical for particles with 
Z = 1 or 2. These distributions have been discussed by Symon ( Symon , 
1948). However, for Z > 2, the distributions are more nearly Gaussian. 
The asymmetry of the distributions arises from statistical fluctuations 
of the relatively small number of collisions in which the energy trans- 
fer is comparable to the average energy loss. The average energy loss 

2 

is proportional to Z , while the maximum energy transferable to an 
electron in a single collision is a function only of the velocity of 
the incident particle: 

« IW > 

hence independent of Z. The fluctuations due to these few high energy 
collisions therefore become less important with increasing charge. In 
addition, variations in detector thickness and particle incidence 
angles, and electronic noise in the amplifiers, tend to broaden the 
distributions. 

A sample of the Dl 1 versus D2‘ distribution for boron and carbon 
having 3=1 is shown in Figure III-l. Since Dl 1 (or D1"0 and D2' are 
nearly identical, they show similar distributions, and can thus be 
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Figure 1 1 1-1 a 

Dl' versus D2‘ for events having P c > 6.4 GV, with 
194 < D3 1 < 203, showing most of boron distribution 
and part of carbon. Most of events not part of main 
distributions are assumed to have suffered detector 
edge effects, as discussed in text. Diagonal lines, 
representing constant Dl" + D2' , are used for 
approximate charge assignment for Figure I I 1-2. 
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Figure IH-lb 

D1 1 versus D2' for events having P c > 6.4 GV, with 
204 < D3 ' $ 214, emphasizing carbon distribution and 
showing some nitrogen and oxygen. Otherwise similar 
to Figure Ill-la. 



Dl CHANNEL 


. . \ . ii. in 

. i . \ . i .1 u 

1. . \ 111 2 

1 . \ l 1. Ill 

\ . . l.\ 22 11 

\ . . . \ 1 2 2 11 1 

\ l 1 . .ll\ 222 1 

\. 1. 1 .1 lK 1 23 112 1 

1 l.\. 121.2.1.. 

\ . 13 11 \ . 

.\ 1 1 . 111131311 IS. . 1 

. \ 1 . 11642 12 2 \2 1 

.1 \l 11 222464426612 \ 

2 . \ i 3. 15446AB6A52 1 ,\ 2 

1 21. 1 \ 1. 1236ACBEAB852 . \ 

1 . 11 \1121 2 119BDBAD66414 1. \ 

I .1 \ .1I218A8BECD88741 . \ 

111 1 \. 1 33AABCCBB67 5 . \ 

...2.1.1....1. l\ 1.38799CAB562... I is 

l 1 .\l22258ECA43642 1 . 

1 .1 2 . \ 144976636 11 

II . I 12 1 1 \ 3342251 1 1 1. 

.1 . 1\ 122122 1 1. 

. 1 1 \ 3111 l . 

s ,1 i . ire 11 1 

\ 1111. 11 \l • 

\ l . 2\. 11 

\. .Ill \ . 

,...L....ii i . .2 . . A. 

\ l 1 . 12 l.l\ 1 


mill 


D2 CHANNEL 

Figure Ill-lb 




30 


combined directly in a number of ways. 'Figure I I 1-2 shows some of the 
distributions of the difference Dl"-D2' for various charges. For this 
purpose, the charges were assigned simply by selecting intervals of 
Dl"+D2' as implied by Figure III-l , and do not represent a rigorous 
charge identification. These distributions are centered at zero (±0.4), 
confirming the validity of Equation II-6. Notice that the central parts 
of the distributions are essentially Gaussian, but there are significant 
numbers of events with Dl"-D2' too large. These may occur through two 
mechanisms: nuclear interactions in the detectors, and edge effects. 

Calculations using geometrical cross sections indicate that <v’1.3% 
(lithium) to ^1.9% (magnesium) of incident nuclei will suffer nuclear 
interactions in D1 ' or D2'. Some of these will be rejected because they 
trigger D4'. The others will have some degree of inconsistency between 
their Dl" and D2‘ pulse heights, which shows up as extra events in the 
wings of the distributions shown in Figure III-2, or as events far from 
the peaks in Figure III-l . 

A much more important source of this background is edge effects 
in DV and D2'. This problem arises from non-uniformity in the electric 
field near the edges of the detectors, causing incomplete collection of 
electron- hole pairs formed by the incident particle. Thus, a particle 
which passes through the edge of Dl' , and then closer to the centers 
of D2 ' and D3‘, will leave a normal signal in D2 ' and D3 1 , but the Dl' 
signal may be anywhere between zero and the normal value. Likewise, a 
particle may leave a normal signal in DV , pass through the edge of D2 ' , 
and then still leave a normal signal in D3', since D3 1 is much larger in 
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Figure III-2 

Distributions of difference DV-D2' for elements 
lithium through carbon. Charge assignments were 
made by selecting contiguous intervals in Dl" + D2', 
and are not intended to be precise. 






diameter than D2'. The geometrical factor for the latter process 
(^0.34 cm sr) is substantially larger than for the former (^.15 cm sr). 
This is borne out by examination of Figure I I 1-2 . The large number of 
events in the wings of the "boron" distribution are due to carbon nuclei 
which passed through the edge of DV or D2‘. This can be seen in 
Figure III-l as events extending towards low D1 ' or D2 ' from the carbon 
peak. It is apparent that there are more events with low D2‘ (high 
Dl"-D2') than vice versa. Since this is purely a geometrical effect, we 
do not affect the elemental abundance ratios by ignoring these events, 
and the absolute fluxes will be correct if the geometry is calculated 
correctly. 

In order to reject the maximum number of background events while 
retaining most of the good events, we need an accurate estimate of the 
shape of the DV-D2' distributions. To obtain this, we selected the 
central portions of the distributions, i.e . 3 |D1"-D2' | < 10, and plotted 
them on probability paper. The results are shown in Figure I I 1-3. These 
plots indicate that the centers of the distributions can be represented 
by Gaussians, and the widths of the Gaussians can be determined fairly 
accurately. Widths read directly from Figure I I 1-3 are shown in 
Figure II 1-4. 

We expect these widths to be composed of a part proportional to Z, 
due to the fluctuations in energy loss, and a part independent of Z, due 
to electronic noise. Thus we should have 
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Figure III-3 

Integrals of Dl"-D2' distributions for |D1"-D2'|<10 
plotted on probability paper for purpose of deter- 
mining widths of distributions. 



FRACTION OF. EVENTS HAVING Dl"-D2'< A 
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A (Channels) 


Figure I I 1-3 
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Figure' II 1-4 

1 

Widths (a) of Gaussian distributions of Dl"-D2' 
obtained from Figure III-3. The dashed line is a 
weighted least-squares fit to Equation (III-2). The 
solid line is the approximation used for event 
selection. 
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A weighted least-squares fit of this function to the observed values is 
shown as the dashed line in Figure III-4. The solid Tine is an 
approximation of the form 

a = o Q + a-|Z (II 1-3 ) 

i 

which appears to fit equally well. The difference between (III-2) and 
( 1 1 1-3 ) is smaller than the uncertainties of the points in Figure 1 1 1-4 , 
and is not significant in any case, as discussed below. 

Figure I I 1-2 indicates that the bulk of the good events have 
| D1 "-D2 ‘ | s 10. Using this as a guide, we chose a tolerance of 2.5a 
in the D1"-D2 I distributions as an event selection criterion. Using 
Figure 1 1-7 to relate Z to Dl l, +D2', and obtaining a from the straight 
line in Figure 1 1 1-4, we require 

|P1"-D2' | < 7.375 + ✓.186 (D1 "+D2' -236) , (III-4) 

i.e. 3 | D1 "-D2 ' j < 2.5a. At this tolerance, we reject only 1.24% 
(nominally) of the events in a Gaussian distribution. Assuming that 
the correct a values are the points plotted in Figure I I 1-4 rather than 
the dashed line, the rejection varies from 0.5% to 1.9%. Therefore, we 
claim that the event selection process is nearly independent of charge, 
and that systematic errors in the elemental abundance ratios from this 
source are at the <1% level. 



Events which are outside the 2.5a tolerance are about 19% of the 
total. As discussed' above, most of these are the result of detector 
edge effects, and can be accounted for geometrically, independent of 
charge. Using geometrical cross-sections, we expect nuclear interactions 
in DV and D2', which have a weak charge dependence, to affect 1% to 2% 
of the total events. Rejection of nearly all of these by (II 1-4 ) , 
introduces a weakly charge-dependent, systematic error at the 1% level. 
Since statistical errors we will be discussing later are ^5%, we feel 
safe in neglecting the above-mentioned errors in the subsequent 
discussion. 

A small fraction of events have similar pulse heights after passing 

through the edges of both DV and D2‘. The geometrical factor for this 
2 

is about 0.16 cm sr. Those falling below the diagonal line in Figure 


1 1-5 are rejected. 
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B. Charge, Velocity Assignment 

Equations (II-2) and (11-5), combined with calibrations derived 
from Figures I 1-3 , I 1-4 , and I 1-7, predict average pulse heights for 
particles with given charge and velocity. An example of this is shown 
in Figure III-5. Here, data taken at cutoff rigidities below 1.8 GV are 
compared with the expected average positions for each charge and velocity 
in an AVG' versus D3 ' matrix (the term AVG' refers to the average of Dl" 
and D2' channel numbers for each event). Events having D3 ' <179 were 
moved to the interval 173 < D3 1 < 179, using the calibration shown in 
Figure I 1-4. This effectively produces a linear response for the entire 
range of pulse heights, at the expense of losing resolution for Z < 3. 

The pair of equations (II-2) and ( 1 1-5) can be solved for charge 
and velocity, given the appropriate pulse heights and calibrations. If 
we allow fractional charges, to account for statistical fluctuations, 
we can assign a unique charge and velocity to each (AVG 1 , D3') pair. 
However, since each channel has a finite width, the following procedure 
was used. Each (AVG 1 , D3‘) pair was considered a rectangle, containing 
a certain number of events. Lines of constant charge were laid across 
this rectangle, and the events were assigned to various charge bins in 
proportion to the area assignable to each bin. Events which, due to 
statistical fluctuations, lie in channels beyond the 6 - 1 line shown 
in Figure 1 1 1-5 , are moved first to the 6-1 line in a direction 
related to the ratio of the widths of the AVG' and D3 ' distributions. 
Events were accumulated in bins 0.1 charge unit wide. Results of 
applying this procedure to the data in Figure II 1-5 are shown in 
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Figure III-5 

Scatter plot of events recorded at P c < 1.8 GV, otherwise 
similar to Figure I 1-5. Calculated average response for 
each charge is indicated by solid lines. Diagonal line 
is fit to peaks shown in Figure II-5. Step-shaped line 
delimits events rejected for AVG' versus D3‘ inconsistency 
as discussed in text. 
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Figure I I 1-6. Similar results for data collected at cutoff rigidities 
above 1.8 GV are shown in Figure III-7. 

These distributions are fit quite well by a sum of Gaussians: 

g 

N(z) = z A. exp t(z - z.) 2 /2a. 2 ] + r(s - z) (III-5) 

3=3 J J J 

for 2.5 s z ^ 9.3. Here z is a continuous variable which represents 
charge. The last term is used to account for the small helium background 
under the lithium peak. The parameter r is a constant for z < s, and 
zero for z > s. In order to allow for possible inaccuracy in detector 
calibrations, z. was allowed to vary quadratically: 

J 

z. = aj + b + gj 2 ( I I 1-6 ) 

J 

Likewise, 

o. = cj + d + hj 2 ( 1 1 1-7 ) 

In all the fits, a, b, c, d, g, h, r, s, and the A. were free parameters. 

J 

As might be expected from looking at Figures 1 1 1-6 and 1 1 1-7 , the results 
were always a = 1 , b= ; c = g= : h = 0, and d = 0.2; i.e. , the charge 
resolution was about 0.2 charge unit, nearly independent of charge. 
Results for each fit are shown in Table III-l. Elements with Z > 9 were 
not fit with this function ( 1 1 1 - 5 ) because D3‘ saturates just beyond 
Z = 9, and the distribution widths are determined only by fluctuations 

l 
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Figure III-6 

Charge histogram indicating assignment of charge to 
events shown in Figure III-5, i.e., events with 
P c < 1.8 GV. Method of assigning charge is discussed 
in text. Fitted curve is Equation I I 1-5. 
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Figure III-7 

Charge histograms for events with P > 1.8 GV. 
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Table III-l. Results of Gaussian Fits (Equation III-5). 


p c> 6V 

0-1.8 

1. 8-4.0 

4. 0-6. 4 

6.4-10.0 

10.0-15.2 

A 3 

106.4 

42.2 

28.6 

24.2 

27.2 

a 4 

63.1 

25.7 

14.6 

12.1 

13.7 

*5 

164.1 

55.3 

40.2 

31.0 

31.6 

*6 

475.4 

189.3 

114.4 

100.5 

118.2 

A ? 

137.4 

53.9 

33.2 

23.0 

29.6 

A 8 

413.6 

163.9 

105.6 

87.5 

118.1 

a 9 

10.4 

3.2 

1.8 

1.2 

3.4 

a 

1.0000 

0.9997 

0.9997 

1.0023 

1.0019 

b 

-0.0266 

-0.0319 

-0.0340 

-0.0683 

-0.0180 

c 

-0.0003 

-0.0006 

-0.0006 

-0.0015 

-0.0010 

d 

0.1987 

0.2020 

0.2032 

0.2182 

0.2252 

10 3 g 

0.378 

0.301 

0.200 

0.956 

0.550 

10 3 h 

0.502 

0.368 

0.366 

0.165 

-0.025 

r 

13.15 

11.53 

8.23 

2.03 

4.62 

s 

3.20 

3.10 

3.10 

3.50 

3.50 

Y 2 

1.882 

1.022 

0.829 

0.794 

0.716 

(x 2 is 

the "reduced 

chi -squared" 

for the fit. 

given by x 2 /v. 

where v=55 


is the number of degrees of freedom.) 
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in AVG'. Likewise, the fits were not extended below Z = 3 because of 
the gain changes of the detectors (See Figures 1 1-3 and 1 1-4 ) , and the 
departure from Gaussian fluctuations for lower charges. 
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C. Abundance Ratios 

For elements 3 < Z < 10, relative abundances can be determined 
with good accuracy from charge histograms such as those shown in 
Figures I I 1-6 and III-7. Such histograms were formed for several ranges 
of cutoff rigidity, and the function defined by Equations ( I I 1-5) 
(III-6), and ( 1 1 1-7 ) was fit to each using an iterative least-squares 
technique. The parameters c, d, and h of Equation (III-7) were also 
determined from a fit to all the data, and the individual data sets 
were fit again using the latter values of c, d, and h. The results of 
both procedures were consistent to within a few events for each charge, 
indicating that variations in the fitted a values from one data set to 
the next are not significant, and have no effect on the relative 
abundances. 

The number of events of each charge may be determined by inte- 
grating the appropriate Gaussian: 

f°° 

n = A. exp[(z - z.) 2 /2o. 2 ] dz = A .a . (III-8) 

J I J • J J J J 

— oo 

and rounding to the nearest integer. We assume that the statistical 
accuracy of this procedure is limited by the Poisson statistics of the 

number of events thus derived. These numbers and key ratios are shown 

/ 

in Table II 1-2. 

Since it was found that this procedure underestimated the area 
under the charge histograms, the following method was used to obtain 


Table II 1-2 . Number of Events Assigned to Each Charge, from Gaussian Fits. 
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final charge assignments. The Gaussian functions derived for all 
charges, and the linear function for helium, were evaluated at each 
charge bin. The events in that bin were then assigned to the various 
charges in proportion to the relative values of the distribution func- 
tions. Thus, essentially all of the events near the center of distri- 
bution were assigned to that charge, while those falling between two 
distributions were split proportionately between them. This method 
assigns all the events in the histograms, and corrects for the possibi- 
lity that Equations III-5, III-6, and III-7 may not be the best repre- 
sentation for the data. The results from this procedure are shown in 
Table I I 1—3 . 

This table also shows helium abundances, which were obtained in a 
different way. For D1 1 and D2' pulse heights below the gain-change point 
(Figure I 1-3) • we used the smaller of the two pulse heights as the 
measurement of the energy loss. This provides a considerable improvement 
in resolution by suppressing the long tail of the Symon distribution 
discussed earlier. For helium, it also has the advantage of removing 
most of the edge-effect events from the distribution. In addition, we 
require |D1'-D2'| < 50 channels, to discriminate against background 
caused by nuclear interactions and extreme Symon fluctuations. About 
2 % of the events are rejected by this method. The helium is quite 
clearly resolved from the singly charged distribution in the resulting 
MIN' versus D3 ' matrix (similar to Figure I 1-5), except at the lowest 
velocities, where the helium distribution overlaps the very large Symon 
tail from the singly charged distribution. For this reason, the helium 


Table 1 1 1-3 . Number of Events Assigned to Each Charge, from Final Algorithm. 
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abundance for P„ < 1.8 GV is uncertain by about 5 - 8% in addition to 
the statistical uncertainty. For P c :> 1.8 GV this presents no problem, 
since the geomagnetic cutoff prevents low-velocity helium from entering 
the telescope. Determination of the helium abundance is therefore 
simply a matter of counting the events within a specified region of 
MIN' versus D3‘ space. 
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IV. Results and Interpretation 
A. Integral Rigidity Spectra 

In order to obtain rigidity spectra, we make use of counting-rate 

v 

information supplied by the Cerenkov rate scaler. This scaler counts 
all events which satisfy the D2'D3 , D4' coincidence requirement, and is 
read out once every 432 msec. During this period, a maximum of three 
events can be analyzed by the pulse-height analysis system. However, 
due to the special design of the rate accumulator ( Althouse et at . , 
1967), the maximum counting rate {'x-lO^ sec"^) is limited only by the 
recovery time of the discriminators. 

Because 0G0-6 is in a polar orbit, it spends a certain fraction of 
its time in regions of high trapped-particle flux. Data taken under 
those conditions are deleted as discussed in Appendix B. For the 
remaining data, the following quantities are accumulated in bins of 
cutoff rigidity (P c ) 0.08 GV wide: 

R = number of Cerenkov scaler readouts, 

C = number of coincidences counted by scaler, 

N = number of Cerenkov events analyzed. 

The total flux of particles in any bin is given by 

J = C/.432RAft (cm 2 *sec*ster) _1 ( IV-1 ) 

2 

where Aft = 3.4 cm *ster is the geometrical factor corresponding to the 
D2 , D3'D4 I geometry, including detector edge regions as discussed 


earl ier. 
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The analyzed events can be assigned charges, as discussed in 
Chapter III, or may be rejected because of Dl"-D2' inconsistency. If 
we designate these as n 2 and n R , respectively, then 

N = n R + s n z . (IV-2) 

Since we have determined that there is no significant charge dependence 
in the criteria for event selection, n R should be composed of events of 
all charges, in proportion to their relative flux. Also, it is reason- 
able to assume that events are analyzed in proportion to their 
respective fluxes, so 

n z /l n z = J z /J , (IV-3) 

or 

J z /J = n z /(N - n R ) . (IV-4) 

However, since most of the rejected events are due to geometrical edge 
effects (as discussed earlier), it is equivalent to correct the geomet- 
rical factor to compensate for the rejection of these events. We define 

G = A/l(zn 7 /N) = ATHl - n„/N) . (IV-5) 

Z L K 
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The flux of a given charge component is then given by 

J z = (C/ .432 RG) (n z /N) = Jn z /(N - n R ) . (IV-6) 

There is some uncertainty in this procedure, because the event selection 
criteria described earlier apply only for 3 < Z < 12. Different criteria 
are applied forZ < 3, resulting in a possible systematic uncertainty of 
^2 - 3% in the absolute fluxes, in addition to the uncertainty in An. 
This problem does not affect abundance ratios for 3 < Z < 10, however, 
since 

J Z /J Z' = n z /n Z' * (IV-7) 

An independent estimate of G/An can be obtained by calculating geomet- 
rical factors directly. We assume G is defined by the gold evaporation 

on the surface-barrier detectors. This was determined to give a 

2 

diameter of 2.34 cm ± 3%, yielding G = 2.75 cm *ster ± 6 %. We estimate 
from the construction of the detectors that the edge region may be 
1 - 1.5 mm wide. This would give an effective diameter of 2.54 - 2.64 cm 
or An = 3.2 - 3.5 cm^’Ster. This results in a ratio G/An * 0.78 - 0.86, 
which is compatible with the value 0.81 obtained from Equation IV-5. 

Integral rigidity spectra for 2 < Z < 10 are shown in Figures IV-1 
to IV-3. Error limits in these figures represent only the statistical 
uncertainties associated with n z for each rigidity bin. The curves 
represent the helium spectrum of Figure IV-1 multiplied by a suitable 
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Figure IV-1 

Integral rigidity spectra of helium, oxygen, lithium, 
and fluorine. The statistical errors in the helium 
fluxes are comparable to the size of the symbols. 
Fluxes are averaged over cutoff rigidity bins having 
edges at multiples of 0.64 GV (0.16 GV for He). Note 
that all the points are statistically independent. 

The curves have the shape of the helium spectrum, with 
normalization determined from fits as discussed in the 
text. 








Figure IV-2 


Integral rigidity spectra of beryllium, boron, and 
carbon. The fitted curves are derived from the 
helium spectrum of Figure IV-1. 
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Figure IV-3 

Integral rigidity spectra of nitrogen and neon. The 
fitted curves are derived from the helium spectrum of 
Figure IV-1. 
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constant in each case. The constants were determined from least-squares 

fits of the spectral points shown in the figures to the helium spectral 

points corresponding to the same range of cutoff rigidities. The results 

of these fits are shown in Table IV-1 . The last column of this table 

2 

contains the probabilities that larger values of x would be obtained if 
the experiment were repeated. It is clear from this table and from 
Figures IV-1 to IV-3 that all the spectra have the same shape, within 
the limits of our statistics. 

In Table IV-2 we compare the abundances relative to helium as 

derived above with those obtained by summing the individual abundances 

from Table 1 1 1-3. We use only the data for P > 1.8 GV for the latter 

c 

calculation, because the data for P < 1.8 GV were not included in the 
spectral fits. The excellent agreement of the results obtained by 
these two methods is further support for the claim that there are no 
statistically significant differences in the shapes of the spectra. 
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Table IV-1. Results of Fits to Helium Spectrum. 



Ratio 

x 2 

V 

x 2 

P(>X 2 ) 

Li/He 


11.8 

12 

0.98 

0.47 

Be/ He 


7.5 

11 

0.68 

0.76 

B/He 


14.2 

13 

1.09 

0.36 

C/He 


20.6 

14 

1.47 

0.12 

N/He 


11.5 

11 

1.05 

0.40 

0/He 

.0251 

9.6 

12 

0.80 

0.65 

F/He 


8.2 

6 

1.37 

0.23 

Ne/He 


10.8 

11 

0.98 

0.46 


The number of degrees of freedom for each fit is v; x 2 is x 2 /v. 
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Table IV-2. Abundances Normalized to Helium. 



From Table II 1-3* 

From Table IV-1 

Li /He 

.00640±. 00025 

.0064 

Be/He 

.00356±. 0001 9 

.0036 

B/He 

.00827±. 00028 

.0081 

C/He 

.0271 +.0005 

.0270 

N/He 

.007531.00027 

.0074 

0/He 

.0252 ±.0005 

.0251 

F/He 

.000501.00007 

.0005 


*Sums for 1.8 < P Q < 15.2 GV 
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B. Abundance Ratios and Interstellar Propagation 
The purpose of making these measurements of the nuclear component 
of the cosmic radiation is to gain information concerning the sources, 
confinement, and propagation of the radiation in the galaxy. At the 
present time, the simplest model which has had reasonable success in fit- 
ting experimental data is the steady-state model of Cowsik et al. 

(1967), in which the assumptions of steady-state source injection, rapid 
diffusion, and slow leakage from the galaxy produce an exponential dis- 
tribution of cosmic-ray path lengths. Inclusion of the effects of ioni- 
zation energy loss, nuclear interactions, and the decay of 10 Be to 10 B 
complete the picture. 

The formulation we use is that of Meneguzzi et al. (1971), which is 
somewhat similar to that of Gloeckler and Jokipii (1969). Following 
Meneguzzi et al. 3 we write the continuity equation: 

+ Ie [b i (E)N 1 (E)D 

- + ‘ , pi (E)n H ]v i N i (E > 

+ + (IV - 8) 

J 1 


N i (E) N i (E) N k (E) 


7+ Q i (E) 
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where we have made the usual assumption that secondaries produced in 
nuclear interactions have the same energy/nucleon as the primary. The 
first term replaces adiffusion term in the "leaky-box" approximation: 
x g is the mean lifetime against escape from the galaxy. The second term 
represents loss of i- type particles by radioactive decay with mean life 
r.j ; the third term, which we have added, represents production of i-type 
nuclei through k -*■ i decay with mean life t^. . Of course, these terms 
are used only for ^Be -* decay in the : range of charges we are con- 
sidering. (E) is the source injection spectrum. The rate of energy 
loss due to ionization is b.. (E) = dE/dt = v^pdE/dx, where p is the 
density of the gas: n^ + n^e^He* The cross sections for destruction 

are a . for destruction on helium and a . for destruction on hydrogen. 

ai pi 

Partial cross sections for production are given by o^. (for j + a -*■ i ) 
and apji (for j + p ■+ i). These include production of unstable nuclei 
which decay rapidly (t«10 y) to i-type nuclei. Finally, v^ is the 
particle velocity; n^ e and n^ are the number densities of helium and 
hydrogen in the interstellar gas. These equations can be solved 
( Gloeoklev and Jokipii, 1969; Meneguzzi et at. , 1971), giving a triangu- 
lar system of equations in the N^. The equation for iron is solved first 
(iron is assumed to be entirely primary), then the equation for each 
lighter nuclide is solved in turn. 

We have made the modification of including as a source of 10 B all 
the 1( ^Be that decays. This is a relatively small correction, amounting 
to about 10% of the flux, or about 3% of the total boron; but it is 
neglected by Meneguzzi et al. In addition, we have used the latest 
value for the ^Be half-life: 1.5 x 10^ y [Yiou and Raisbeck, 1972). 
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The model includes nuclear interactions with interstellar helium as 

well as hydrogen. We assume that the interstellar gas is composed of one 

hydrogen atom and 0.1 helium atom per cubic centimeter. The assumed 

source composition is shown in Table IV-3. This is basically the same as 

4 14 

that used by Meneguzzi et al. t but we have adjusted the He, N, and 
16 0 abundances slightly to fit our data. The source spectra are assumed 
to be power laws in total energy, with exponent -2.6. The escape mean 
free path, A , is a free parameter. > 

The result of the calculation is a differential energy spectrum for 
each isotope from to ^0, plus 28 Ne, 2 Slg, 28 Si , and 56 Fe. To com- 
pare these results with our data, we convert each differential energy 
spectrum to an integral rigidity spectrum, sum the resulting fluxes over 
all the isotopes of each element, and divide by the oxygen flux. The 
results are shown as the dashed lines in Figures IV-4 and IV-5; the data 
points are taken from Table III-3. Abundance ratios shown are calculated 
for A g = 6, 8, and 10 g/cm . Neglecting interstellar helium would have 

resulted in higher predicted ratios, or lower values of a q deduced from 

2 

the data. The difference in A Q is approximately 1 g/cm . The curves 

2 

shown in the figures would thus correspond to A g = 5, 7, and 9 g/cm if 
the helium were not included in the calculations. Our calculated L/M for 
A g = 6 g/cm would agree very closely with that calculated by Mason 
(1972) if we neglected the interstellar helium. 

We cannot compare the data directly to these curves, because the 
modulation of galactic cosmic rays by the solar wind may be expected to 
alter the composition. This problem is particularly severe for our data 
which were obtained during solar maximum. The usual procedure for 
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Table IV-3. Assumed Cosmic Ray Source Abundances. 


Isotope 

Abundance (this work) 

SST 

WDK 

4 He 

32±1 

27±5 

- 

12 c 

l f 

1 + 

1 + 

14 n 

0.07* 

0.12±.03 

0.110±.016 

16 0 

1.2±0.1 

1.02±.06 

1.06 

20 ne 

0.25** 

0.20±.03 

0.17 

24 Mg 

0.29** 

0.27±.04 

0.23 

28 Si 

0.26** 

0.23±.04 

0.17 

56 Fe 

0.25** 

0.23±.05 

0.19 


+ Normal ization 

*This value was used to calculate the curves in Figure IV-5. 
The data appear to imply a value of 0. 10+.01 for N/C at the 
source. 

**Abundances of elements heavier than oxygen are taken from 
Meneguzzi et al. (1971), who adjusted them to compensate for 
the neglect of secondary 1 4<Z<26 . 

SST - Shapiro et al. (1971) 

WDK - Webber et al. (1972) 



Figure IV-4 


Ratios Li/0, Be/0, and B/0 with predicted interstell 

values (dashed lines), and values at 1 A.U. (solid 

lines). Numbers at right are values of escape mean 

2 

free path, A , in g/cm . Neglecting interstellar 

helium yields similar curves with A g decreased by 
2 

~1 g/cm for each curve. 
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Figure IV-4 
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Figure IV-5 

Ratios He/0, C/0, and N/0. Source abundances were 
adjusted to fit data. 



NITROGEN/OXYGEN CARBON/OXYGEN HELIUM/OXYGEN 
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Figure IV-5 
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attacking this problem is to make some assumptions about the functional 
form of the modulation, estimate its effect on the abundance ratios, and 
extrapolate the observations to interstellar space. 

We take the opposite approach here: the fluxes calculated from the 

interstellar propagation model are used as boundary conditions for numer- 
ical solutions of the transport equation governing the motion of cosmic 
rays in the interplanetary medium: 

3U/3t = - *.(VJJ) + i t(^*V ul )3(aTU)/3T] + V*(k-^U) (IV-9) 

W 3 W 

where U is the number of particles per unit volume per unit interval of 

kinetic energy, V is the solar wind velocity, T is the kinetic energy, 
w 

is the diffusion tensor, and a = (T + 2M)/(T + M), where M is the par- 
ticle's rest energy. This description of solar modulation has recently 
been reviewed by Jokipii (1971). 

Equation IV-9 may be simplified by making the following assumptions: 

1) spherical symmetry, 

2) constant, radial solar wind velocity, 

3) no time dependence, 

4) no modulation beyond a certain radial distance, D. 

Using these assumptions, we get ( Garrard , 1973): 

t 

r" 2 V w a(r 2 U)/3r - (2V w /3r)a(aTU)/9T - r" 2 9(r 2 K 9U/9r)/9r =0 (IV-10) 

to be solved with the boundary condition U(T,r=D) = U(T,interstellar). 
Here k is the radial part of the diffusion tensor. 
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The solution we use is that of Garrard (1973), with the following 

parameters: D = 6.1 AU, V w = 400 km/sec, and k = k Q gf(R), where f(R) = R 

for R * R Q , f (R) = /RR q for R < R Q , R q = 1.5 G V, and-k =1.1 x 10 21 
2 - 1-1 

cm *sec *GV . These parameters were determined ( Cummings , 1973; 
Garrard, 1973) for 1969, when most of our data were taken. The transport 
equation is solved separately for each isotope; the resulting fluxes at 
1 A.U. are converted to integral rigidity spectra; these are summed over 
all isotopes of each element and divided by the oxygen spectrum. The 
resulting ratios are shown as the solid lines in Figures IV-4 and IV-5. 

A number of conclusions can be drawn from Figure IV-4, which shows 
the abundances, relative to oxygen, of the "light 2 * * * * * * * * 11 elements, which are 
assumed to be absent from the sources. First, it appears that the 

2 

beryllium/oxygen and boron/oxygen ratios both correspond to A g - 10 g/cm 

2 

at ^2 GV, decreasing to 6 or 7 g/cm at 10-15 GV. (The boron point at 

1.8-4 GV apparently is low due to a statistical fluctuation, as may be 
seen by examining the boron spectrum in Figure IV-2.) A dependence of 
escape mean free path on rigidity was suggested by Cowsik et al. (1967), 

but has since been neglected by most authors ( e.g., Shapiro et al. , 1970; 

Webber et al. , 1972; Mason, 1972), usually in favor of a model employing 
an exponential distribution with the short path lengths suppressed. 

Von Rosenvinge et al. (1969) conclude that the lack of variation of L/M 

(the commonly used ratio [Li + Be + B]/[C + N + 0]) with energy is 

evidence for a path length distribution which is either a delta-function 
or a Gaussian with a very small width. More recently it has been pointed 

out {Mason, 1972) that the constancy of this ratio at low energies 

(s300 MeV/nucleon) is evidence for adiabatic deceleration of the cosmic 
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rays in the expanding solar wind. At intermediate energies, such as 

those we observe, small solar modulation effects may combine with a 

variation in a 0 with rigidity (or energy) to produce a nearly energy- 

independent L/M. At higher energies (;>4 GeV/nucleon) a decrease of 

secondaries/primaries with increasing energy ( Juliusson et at. , 1972; 

Webber et at., 1973; Smith et at., 1973) points strongly towards an 

energy-dependent leakage mean free path. 

The second observation we make from Figure IV-4 is that there 

appears to be too much lithium, if we assume that the A g implied by 

beryllium and boron is correct. There are four possible explanations 

for this: 1) currently assumed cross sections for producing lithium 

assumed 

from heavier cosmic rays are too small, 2) the/ total inelastic cross 

7 

sections for destroying Li and Li are too large, 3) there is some 
thermonuclear 7 Li present in the cosmic ray sources, and 4) our measure- 
ment is wrong. 

The last possibility would seem to be the simplest explanation, 
since it is rather difficult to separate lithium completely from spill- 
over from the very large helium peak. We have examined the helium dis- 
tributions in all our detectors, and carefully corrected the lithium for 
contamination from this source. The probability that 30% of our remain- 
ing lithium is still mis identified helium is rather small, but probably 
non-zero. 

Possibilities (1) and (2) above certainly cannot be ruled out. 
Measurements of fragmentation of L nuclei in nuclear emulsions 
( Cleghom , 1967) indicate that the isotopes of lithium and beryllium may 
be exceptional cases which do not agree with predictions of cross 
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sections based on geometrical or semi -empirical ( e.g . 3 Silberberg and 
Tsao, 1971) formulas. Extensive cross section measurements will be 
required to resolve this point. 

The third possibility is the most interesting from an astrophysical 
point of view. The existence of a source of interstellar 7 Li was pro- 
posed by Meneguzzi et at. (1971) to explain the large 7 Li/^Li ratio 
observed on the earth and moon, and in meteorites and some stars. The 
existence of a substantial amount of 7 Li in cosmic-ray sources would 
seem unlikely, however, because the deduced source abundances of the 
heavier elements indicate that the sources are highly evolved objects, 
such as supernovae, which would probably not contain sufficient lithium 
to account for this observation. Similar discrepancies in the measured 
lithium abundances have appeared previously ( Von Rosenvinge et at. , 

1969; Mason, 1972), indicating either the difficulty of making an 
unambiguous determination of the lithium abundance, or the necessity 
of revising the model. 

In Figure IV-5 we show ratios of elements which presumably are 
present at the cosmic ray sources. As mentioned previously, the source 
abundances have been adjusted to fit the data, although it appears from 
the figure that we may have underestimated the nitrogen abundance in 
the source by about .03. The important feature to be noted here is that 
the helium and nitrogen (and possibly the carbon) also show the tendency 
toward shorter path lengths at higher rigidities. This is independent 
of the assumed source abundances, and requires only that all the source 
spectra have the same shape. 
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2 

We conclude (1) that an escape mean free path varying from MO g/cm 
at *2 GV to -v6 g/cm 2 at MO GV enables the simple "leaky box" model to 
fit the data quite well, and that no modifications to the exponential 
path-length distribution are necessary; and (2) that individual element 
ratios, such as Be/0, are superior to combined ratios, such as L/M, for 
testing propagation models. Of course, individual isotope ratios will 
be even more enlightening when they become available. 

For comparison we show some published results of other authors, with 

our calculations and data for some of the 1 combined ratios. Figure IV-6 

shows the L/M ratio obtained by Mason (1972) on the IMP-5 spacecraft. 

The dotted curve is his calculation based on an exponential path-length 

2 

distribution with a mean of 6 g/cm , which would correspond to A g - 
2 

7 g/cm if helium is included. Some balloon results are shown in 
Figure IV-7. These have been extrapolated to the top of the atmosphere 
by Webber et at. (1973). These results are strong evidence for a con- 
tinuing decrease in the A 0 with increasing energy. 

Our underestimation of the amount of nitrogen in the cosmic ray 
sources, and the ambiguity in the comparison of calculated and measured 
lithium, are examples of the pitfalls inherent in the use of ratios of 
this type. The improvement in statistics obtained by combining all the 
L elements is more than offset by the increased difficulty of inter- 
preting the results correctly. 

Finally, we turn to the problem of estimating the confinement time, 
or "age" of the cosmic rays by examining the effect of ^Be decay. This 
problem has never been completely resolved, because of the poor statis- 
tics involved in measuring beryllium, the uncertainties in the cross 
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Figure IV- 6 

L/M from this experiment (upper panel) and from 

Mason (1972, lower panel). Solid curves are 

predicted values at 1 A.U. , for values of A g 

specified at right. Dotted curve is Mason's 

2 

calculation for A g = 6 g/cm , with no inter- 
stellar helium. 



(Li + Be + B)/(C + N + 0) 
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Figure IV-6 
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Figure IV-7 

L/(C + 0) from this experiment (upper panel) and 
from balloon measurements as indicated. Balloon 
data were extrapolated to the top of the atmosphere 
by Webber et al., who also converted the Berkeley 
data from rigidity to energy intervals. 
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Figure IV- 7 
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sections, and smallness of the expected effect. Some authors have 

claimed or implied substantial decay of ^°Be ( Lezniak et al. 3 1970; 

Webber et al . 3 1972) while others have claimed survival [Von Rosenvinge 

et at. 3 .1969; O'Dell et at., 1971). 

Our results are shown in Figure IV-8, where we have used a value of 
2 

8 g/cm for A g , and varied the average density of the interstellar 

_3 

hydrogen from 0.1 to 10 cm , keeping He/H =0.1. The average confine- 
ment time, t , can be related to A g by 

T e = V Bc <"h H H + VW (IV ’ 11) 

where n and M are the appropriate number densities and masses. For 
n^ = 1 cm -3 , - 1.14 x 10 14 /6 sec * 3.5 x 10 6 /p years. Our data 

appear to favor a value of n H not too different from 1 cm , so a life- 
time of a few million years is implied. Our ratio for Be/B for all data 
with P c >1.8 GV is 0.43 ± 0.03. This is only 0.15 standard deviation 
from the value calculated by O'Dell et al. (1971) for 10 Be survival, but 
is 2.5 standard deviations from their value assuming decay of 10 Be. 

From the Be/B ratio for P c < 1.8 GV, we obtain a mean confinement time 
6 5+0 5 

of 10 ‘ * years. The points at higher rigidities, when averaged 

together, imply 10 7 years. 

There may be some additional effect due to the variation of A g with 

rigidity, but it will not be possible to investigate this until much 

4 

better data are available (at least 10 Be events, preferably observed 
at solar minimum) and until all the appropriate cross sections are 
accurately measured or reliably estimated. 
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Figure IV-8 

Ratio Be/B. Solid curves are predicted values at 

1 A. U. Numbers at right are average density of 

3 

interstellar hydrogen, in a toms /cm . Helium density 

is one-tenth of hydrogen density in each case. 

2 

Curves are calculated using A g = 8 g/cm . 
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V. Summary and Conclusions 

We have made the first comprehensive satellite measurement of the 
abundances and spectra of cosmic ray nuclei using the geomagnetic field 
as a spectral analyzer. This was made possible by the low-altitude 
polar orbit of the 0G0-6 spacecraft,, and its attitude control system 
which kept the Caltech experiment aimed radially away from the earth. 
The advantages of this type of measurement over the more common method 
of using high-altitude balloons are (1 ) the ability to sample a wide 
range of cutoff rigidities, and (2) the elimination of the difficult 
problem of correcting for atmospheric secondaries. The instrument used 
the technique of making simultaneous measurements of energy loss and 

V 

Cerenkov radiation. 

By making a very -careful and extensive study of the properties of 
the instrument, we were able to obtain data with very low background, 
good charge resolution (a = 0.2 charge unit), and negligible charge- 
dependent bias. This was done using a relatively simple instrument, 
containing a minimum of mass which might induce nuclear interactions. 

We found that all the elements 2 < Z < 10 have similar integral 
spectra over the range of cutoff rigidities from ^2 to ^15 GV, 
approaching a p -1 * 6 power law at rigidities *8 GV. Ratios of various 
elements are essentially independent of rigidity over the range 
observed. The values of the ratios and their dependence (or lack 
thereof) on rigidity, combined with calculations of interstellar pro- 
pagation and solar modulation imply the following conclusions: 

1) A simple "leaky box" model for galactic confinement of cosmic 
rays is adequate to explain the data. 
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2) The leakage, expressed as a mean free path, must vary with 

p 

rigidity, having a value decreasing from VI 0 g/cm at V? GV to 
Vi or 7 g/cm^ at V10 GV. 

3) The effect of interstellar helium on the propagation of cosmic 
rays cannot be neglected. 

4) The effect of solar modulation on elemental abundance ratios is 
important for rigidities less than 3 or 4 GV (kinetic energies 
<1 or 2 GeV/nucleon) , and is probably dominant below 200 or 300 
MeV/nucleon, at solar maximum. 

5) Individual element ratios, such as beryll ium/oxygen, are pref- 

erable to combined ratios, such as L/M, for testing propagation 
models. 

6) Accurate measurements of production and destruction cross sec- 
tions for the isotopes of lithium, beryllium, and boron will be 
needed before finer details of the propagation can be studied. 

7) The ratio of beryl! ium/boron implies an average density of 
neutral hydrogen of about 1 /cm and a mean confinement time 
of iq 6 * 5±0 * 5 years. Survival of most of the 10 Be is favored. 
This implies confinement of the cosmic rays mainly within the 
galactic disc. 
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APPENDIX A. Data Processing and Assignment of Cutoff Rigidities 

It is beyond the scope of this thesis to discuss in detail all the 
processing which was done on the data produced by the experiment, since 
most of it was of a bookkeeping nature and had no direct bearing on the 
physics involved. We simply present here an outline of what was done, 
with emphasis on those aspects which have some relevance to the inter- 
pretation of the results. More details may be found in a number of 
reports on the subject {Brown, 1971a, 1971b, 1972). 

The first stages of processing consist of merging experiment data 
and satellite orbit data, checking for and deleting data of poor qual- 
ity, and extracting and editing information relevant to the Cerenkov 
Telescope. This is done primarily to reduce the expense and incon- 
venience of handling large numbers of tapes, and to put the data in a 
more convenient form. 

The next step is to assign vertical cutoff rigidities. Satellite 
orbit data are available every 60 seconds; these are interpolated to 
15-second intervals, then used to calculate cutoff rigidity from a grid 
calculated by Shea and Smart ( Shea et at. , 1965, 1968; Shea and Smart , 
1967; Shea, personal communication, 1971). This grid represents 
"effective" cutoffs, which are derived by averaging over allowed 
rigidities in the penumbra! region ( Shea et at . , 1965). The grid we 
used was modified from a 2° by 2° grid kindly supplied by M. A. Shea 
(personal communication, 1971). The supplied grid was interpolated 
from a 5° by 15° grid previously published ( shea et al., 1968). We 
found it necessary to refine the interpolation in the region bounded by 
latitudes 0°-30°N and longitudes 245°-325°E, using the 5° by 5° grid 
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calculated by Shea and Smart in that region, as well as points calcu- 
lated at one-degree latitude intervals for longitudes 260°E and 285°E. 
This was done by hand, using an iterative graphical technique. 

Since this grid is calculated for an altitude of 20 km, the cutoff 

is extrapolated to spacecraft altitude using the relation {smart and 
-2 

Shea , 1967) P^L where L is the Mcllwain parameter ( Mcltwain , 1961). 
Finally, the cutoff rigidities are interpolated to 144 msec intervals 

• V 

and assigned to each analyzed event and each readout of the Cerenkov 
scaler. 

It is difficult to estimate the accuracy of these calculated ver- 
tical cutoffs, since few direct measurements exist. One recent pair of 
measurements near Palestine, Texas {Pennypacker et al., 1973) using a 
balloon-borne superconducting magnetic spectrometer gives results which 
agree to within 0.1 GV with the corresponding interpolated grid values. 
Since the quoted uncertainty in the measurements is 0.1 GV, we consider 
the agreement to be excellent. However, it remains an open question 
whether discrepancies might be found at other geographic locations, 
although the calculations for high latitudes (P c s 1 GV) are probably 
not reliable ( Fanselow and Stone, 1972) due to the neglect of the field 
due to external currents ( Smart et al., 1969). Also, the calculated 
cutoffs correspond to a quiet geomagnetic field model. We have made no 
attempt to correct for geomagnetic activity, the effect of which is to 
lower the cutoffs from the quiet-time values. This effect may be 
expected to be most pronounced at high latitudes (low cutoffs) resulting 
in an apparent distortion of the integral rigidity spectra. However, 
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because the various abundance ratios are nearly independent of rigidity 
(Figures IV-4, IV-5, and IV-6), distortion of these ratios by fluctua- 
ting cutoffs is a second-order effect, and can be neglected without 
affecting the conclusions of this thesis. 
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APPENDIX B. Identification of Geomagnetically Trapped Particles 
The Caltech Solar and Galactic Cosmic Ray Experiment on 0G0-6 
consists of three separate charged-particle telescopes which share 
pulse height analyzers and readout circuitry. In order to understand 

V 

the response of the Cerenkov Telescope/electronics system, it is neces- 
sary to consider the effects of the other two telescopes. For this 
purpose, it is sufficient to consider the following: 1) the energy 

ranges and geometrical factors of the telescopes, 2) the way in which 
events are selected for analysis, 3) the effects of analyzer dead time, 
4) the characteristics of the particle fluxes seen by the instrument. 
The Range Telescope analyzes electrons with E ;> 1 MeV and nuclei 

with E :> 1 MeV/nucleon, with a geometrical factor decreasing from 
2 2 

M.5 cm -ster at low energies to 0.2 cm «ster at E = 300 MeV/nucleon. 

V 

It has the same priority for event analysis as the Cerenkov Telescope. 
The Flare Telescope analyzes protons (17 to 100 MeV) and alphas (70 to 
400 MeV), with a geometrical factor of about 0.02 cm *ster ( Althouse 
et al. j 1967). The Flare telescope has the highest priority for event 
analysis. The aluminized Mylar window of the Range Telescope failed in 
January, 1970, causing the front detector to respond to sunlight. This 
detector was disabled by ground command, and the logic required for 
event analysis was modified, causing the Range Telescope to analyze 
electrons with E ;> 240 keV after this time. 

The pulse-height analysis system is capable of reading out one 
event every 144 msec. When the spacecraft is in a region of geomag- 
netically trapped particles, or in the polar regions during a solar 
flare, the counting rates of some or all of the telescopes are much 
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higher than this. If the Flare Telescope rate is »7/sec, all the 
analyzed events will be Flare events. If the Flare Telescope rate is 

V 

$7/sec, the Range and Cerenkov Telescopes will compete for whatever 
analyzer time remains according to the ratio of their counting rates. 
Because the Range and Flare Telescopes have lower energy thresholds 
than the Cerenkov Telescope, they usually have much higher counting 
rates in regions of high trapped-particle flux. As a result of this, 
plus the event-analysis priority given to the Flare Telescope, rela- 
tively few Cerenkov Telescope events are analyzed in the trapped- 
particle regions. We take advantage of this fact to identify and 
eliminate data which are contaminated by trapped particles, or by solar 
flare particles in the polar regions. 

The quantities R, C, and N, defined in Chapter IV (number of 
Cerenkov scaler readouts, number of coincidences, number of analyzed 
events) are accumulated in one-minute intervals. Typical values are 
R < 139 min" 1 ; C = 10-40 min” 1 for galactic cosmic rays, and C » 

40 min” 1 for trapped particles; N * .9C for intervals with no trapped 
particles, and N - 0 for intervals containing trapped particles. The 
ratio N/C is a sensitive indicator of the presence of trapped particles, 
and can be used in a "yes-no" fashion to determine whether a one- 
minute interval should be considered contaminated with a significant 
number of trapped particles. For this purpose, any interval having a 
low efficiency for event analysis (N/C < r) was designated as contam- 
inated, where r is a parameter discussed below. 
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The sum of 15 months of data is shown in the following figures. 
Figure B-la shows the behavior of the quantity C/R as a function of cut- 
off rigidity (P c ). Notice the pronounced increases in counting rate at 

P =1.2 GV (outer trapping zone, mostly electrons) and at P_ =9.9 GV 
c c 

(inner trapping zone^ mostly protons). The same quantity is shown in 
Figure B-lb, where a distinction has been made between readouts having 
N/C < r and those having N/C > r. These are designated as "low 
analysis-efficiency" data (L) and "high analysis-efficiency" data (H) 
respectively. 

Likewise, the ratio N/C is shown in Figure B-2a for all the data, 
and in Figure B- 2b after making the separation. Notice that N/C is 
practically constant in the "high efficiency" data. The small varia- 
tions are due to the use of analyzer time by the other telescopes, and 
are not directly correlated with the Serenkov rate. A value of 0.75 
for the parameter r was chosen, after some trial and error, to provide 
the best separation of the behavior of high efficiency and low effi- 
ciency C/R and N/C. Lower values allow too much contamination in the 
high efficiency data; and higher values cause the rejection of too 
much good data. However, any value of r Z 0.25 rejects the vast 
majority of contaminated intervals. 

Because the maximum counting rate for galactic fluxes is about 
40 min"\ and the maximum number of analyzed events is about 417 min \ 
the probability is less than 10% that a Cerenkov event will not be 
analyzed because the analyzers are busy with a previous Cerenkov event. 
Thus, the probability of failing to analyze ten (25%) of the Cerenkov 
events is less than 0.1^ = 10"^. However, since some of the analyzer 
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Figure B-1 

Number of D2'D3'M' coincidences per scaler readout 
observed as a function of cutoff rigidity. Rigidity 
bins are 0.08 GV wide. 

a. Total of all observations; 

b. Observations having high efficiency for 
event analysis (H) separated from low 
efficiency data (L). Large fluctuations 
in the (L) data for P c z 13 GV are due to 
the small number of L-type readouts at 
high cutoffs. 
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Figure B-2 

Fraction of Cerenkov events analyzed as a function 
of cutoff rigidity. 

a. Total of all data; 

b. Low efficiency (L) and high efficiency (H) 
data shown separately. 
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capacity is being used by the other telescopes, the actual probability 
is much higher, resulting in the loss of about 10% of the Cerenkov 
events, as indicated in the H-curve of Figure B-2b. However, this loss 
rate is independent of C/R for the high efficiency data. 

It can be seen in Figure B-lb that there is some contamination of 
the high efficiency data by trapped particles at cutoff rigidities 
between 11 and 12 GV. This must be composed entirely of singly charged 
particles, since no similar "bump" appears in the helium spectrum of 
Figure IV-1. In fact, the shape of the total "spectrum" implied by the 
H-curve in Figure B-lb is substantially different from that of the 
helium spectrum, due to the large flux of return albedo at high (*2 GV) 
cutoffs. For this reason we cannot obtain a spectrum of primary pro- 
tons. However, the spectra for all Z > 1 particles are derived cor- 
rectly using the method described in Chapter IV, because the albedo 
and trapped particle fluxes are singly charged. Splash albedo has been 
measured directly ( Wenzel , 1968) at Churchill, Canada (P c « 1 GV) and 
Palestine, Texas (P =4.5 GV) using an instrument essentially identical 
to that on OGO-6. The measurements of helium in the splash albedo 

(T. Garrard, personal communication, 1973) imply an upper limit of 
-2 -1 -1 

sl m *sec *ster for the integral flux above 300 MeV/nucleon, assum- 

_2 

ing a differential energy spectrum «E . The spectrum is very likely 

\ 

steeper than that above 300 MeV/nucleon, implying an even smaller upper 
limit. These considerations indicate that the flux of albedo helium 
should be «0.6% of the primary helium flux, causing no ambiguity in 
the derivation of the primary spectrum. The same considerations pre- 
sumably apply for heavier nuclei as well. An independent verification 
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of this conclusion can be obtained by examining event distributions at 
high cutoffs ( e.g Figure I 1-5 ) and noticing that there are no events 
with low velocities. 

The number of scaler readouts R is shown in Figures B-3a, b, and c. 
The high efficiency readouts in Figure B-3b represent the sampling time 
available for the analysis done in this thesis. 
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Figure B-3 

Observing time as a function of cutoff rigidity. 

a. Total of all observations; 

b. Readouts having high efficiency for event 
analysis; 

c. Readouts having low efficiency for event 
analysis. 
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